Superconductivity near lattice instability: the case of NbCi-^N^ and NbN 
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Using density functional theory within the local density approximation we report a study of 
the electron-phonon coupling in NbCi-^N^ and NbN crystals in the rocksalt structure. The Fermi 
surface of these systems exhibits important nesting. The associated Kohn anomaly greatly increases 
the electron-phonon coupling and induces a structural instability when the electronic density of 
states reaches a critical value. Our results reproduce the observed rise in T c from 11.2 K to 17.3 K 
as the nitrogen doping is increased in NbCi-zN^. To further understand the contribution of the 
structural instability to the rise of the superconducting temperature, we develop a model for the 
Eliashberg spectral function in which the effect of the unstable phonons is set apart. We show that 
this model together with the McMillan formula can reproduce the increase of T c near the structural 
phase transition. 



I. INTRODUCTION 

The discovery of superconductivity in boron-doped 
diamond renewed interest in phonon-mcdiated type 
superconductors. 1 Some of these materials like CS3C60 2 
and MgB 2 3 are reported with transition temperature T c 
up to almost 40 K. In the case of MgB2, it has been 
suggested that the Kohn anomaly is responsible for the 
enhanced electron-phonon coupling leading to a larger 
critical temperature T^This phenomenon has been ar- 
gued to raise electron-phonon coupling in copper oxide 
superconductorsP In the case of transition metal car- 
bides, the relatively high T c is also explained by a Kohn 
anomalyP The nitride NbN in the rocksalt structure is 
commonly reported as having the highest T c among the 
carbides and nitrides at 17.3 KPEl From a technological 
point of view, this material is interesting because of its 
possible application to induce superconductivity in car- 
bon nanotube junctions.^ Recent phonon calculations^ 
showed that the rocksalt phase of NbN is unstable. It is 
known from experiment, however, that the rocksalt phase 
can be stabilized in the alloy NbCi-^N^ and in the nitro- 
gen deficient NbN crystal!^ In the present paper, we re- 
port ab initio calculations of electron-phonon coupling in 
NbCi-^Nz and in the nitrogen deficient NbN to establish 
a parallel between the Kohn anomaly and T c . We elabo- 
rate a model based on the Eliashberg spectral function, 
a 2 F(ui), explaining the enhancement of electron-phonon 
coupling due to the Kohn anomaly near the structure 
phase transition. 



II. METHOD 

The calculations reported in the present study were 
carried out using density functional theory (DFT) 
within the local density approximation (LDA) as im- 
plemented in the ABINIT codeP^ Norm- conserving 
Trouiller-Martins pseudopotentials were used to repre- 



sent the interaction of the valence electrons with the 
atomic cores. The sampling of the Brillouin zone was 
done using a 24x24x24 fc-grid with a gaussian broadening 
of the Fermi-Dirac distribution of 5 mHa. Wave functions 
of the electrons are expanded with a plane-wave basis up 
to an energy of 35 Ha. 

The NbCi-^Na; crystal was calculated within the 
virtual-crystal approximation (VCA) by simple mixing 
of the ioni c pseudopotentials of carbon and nitrogen 
atomsPSEH The nitrogen deficient NbN crystal was sim- 
ulated by removing x electrons per unit cell from a com- 
plete NbN crystal which we will denote by NbN 31 . A 
uniform background charge maintains charge neutrality. 
Phonon spectra and electron-phonon coupling were eval- 
uated within a linear response theorjttSLLS] on a 12x12x12 
q-points sublattice. 

III. RESULTS 

The calculated density of states (DOS) is reported in 
Table [T] In another study, Isaev et al. published results 
of ab initio calculations for NbCP^The main differences 
between their study and the present one are their use of 
ultrasoft pseudopotentials with the generalized gradient 
approximation (GGA) PBE for the exchange-correlation 
functional. We also used a denser grid of fc-points for 
integration in the Brillouin zone (24x24x24 compared to 
18x18x18) and a denser g-points grid for the phonon cal- 
culations (12x12x12 compared to 8x8x8). Our results 
agree with theirs to the expected precision of this type 
of approach given the different parameters used. 

The experimental data reported by TotlP are in agree- 
ment with our results given in Table [i] In particular, the 
experimental value of the DOS at the Fermi level, Np, 
for NbC is about 0.2 eV _1 per atom, which corresponds 
to 0.4 eV _1 per unit cell, which is within 10% of our cal- 
culated value of 0.361 eV _1 per unit cell. This agreement 
is satisfactory considering that the experimental value is 
obtained through specific heat measurements using the 



2 




FIG. 1. (color online). Fermi surface for NbC rendered by 
XCrysDenP^ The cylinders exhibit important nesting for k in 
the TX direction. Substituting carbon by nitrogen increases 
the radius of the cylinders. 



Sommerfeld free electron model. Furthermore, the de- 
duced A from experimental data quoted by Toth employs 
the version of the McMilllan formula using the Debye 
temperature which is extracted again from the specific 
heat. Hence, a direct quantitative comparison with our 
results is not possible but we note that Toth quotes a 
value of 0.6 for A whereas we obtain 0.682. Considering 
the differences mentioned above, this correspondence is 
acceptable and validates our approach. 

Upon inspection of the NbC Fermi surface, we see that 
it is made of arms lying along the TX axis as depicted in 
Fig. [TJ This topology, which allows important nesting for 
q vectors connecting opposite faces of the arms, results 
in Kohn anomalies. This is similar to the case of TaC re- 
ported by Noff singer et alPFor NbCi-ajN^, adding elec- 
trons by partially substituting C by N atoms increases the 
radius of the arms of the Fermi surface which leads to en- 
hanced nesting due to a larger phase space. This change 
of the Fermi surface shifts the wavevector of the resulting 
softened phonon towards larger q. This phenomenon is 
readily observed in the phonon spectrum along TX di- 
rections as shown in Fig. [3] for the longitudinal acoustic 
(LA) branch. 



The softening observed in the NbC spectrum, Fig. 2(a) 



is consistent with our analysis of the Fermi surface men- 
tioned above. For NbCo.sNo.s and NbN - 4 , the softening 
becomes an instability with imaginary phonon frequen- 
cies as illustrated in Fig. |2(c)' and 2(d) respectively. We 



see from Fig. 2(b) that NbN in the rocksalt structure is 



indeed unstable. We will refer to the density of states at 
the Fermi level of the structure at which the softening be- 
comes an instability as N c . The calculated values for this 
quantity arc expected to be lower than the experimental 
values^ because we do not take into account anharmonic 
effects which are known to diminish the softening.^ The 
quantity N c will prove to be fundamental in the elabo- 
ration of our theoretical model for the electron-phonon 
coupling. 

Table [TJ reports the values, calculated within linear re- 



100 



so 



> 

a 6o 



40 



(a) NbC 




100 



so 



> 

a 60 



40 




20 



100 



80 



(b) NbN 





> 

a 6o 



40 



(c)NbC _ 6 N a4 




100 



so 



> 

e 60 



40 



(d) NbN 




FIG. 2. Phonon band structure for (a) NbC, (b) NbN, (c) 
NbCo. 6 No.4, (d) unstable NbN - 4 
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TABLE I. Density of states at the Fermi level, Nf, and results of linear response theory calculations for the electron-phonon 
coupling in NbCi-^Nz and in nitrogen deficient NbN. The total electron-phonon coupling constant, A, and the contribution of 
the acoustic branches, A ac , are reported separately. The weighted average of the phonon frequency, cJi og , is needed to compute 
the transition temperature T c using the McMillan formula |T]). The parameter fi* is set to 0.1 in order to fit the NbC T c to its 
experimental value of 11.2 Kp^The experimental values of T c are also reported for comparison.^ 
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FIG. 3. (color online). Phonon dispersion along the TX di- 
rection for various concentration on N. The evolution of the 
Kohn anomaly is clearly seen. The softening becomes more 
important with the number of electrons and is also shifted 
towards X. 



tween calculations and experiments^ but this has little 
effect on the model elaborated below. For NbCi-^Nj., 
the rise of T c with the increase of the nitrogen concen- 
tration obtained in our calculations reproduces the trend 
in the experimental data^l The increase in the theoret- 
ical values is larger than the experimental one but the 
agreement is satisfactory considering the approximations 
involved in the McMillan equation. For NbN, the only 
experimental T c reported is 17.3 K.^l According to our 
study, the NbN in the rocksalt phase is only stable when 
the system is charged. We suspect that the experimen- 
tal system is non-stoichiometric with deficient nitrogen. 
Therefore, it is appropriate to compare the experimental 
value with our NbN calculations near the structural in- 
stability. We see from Table [I] that we obtained a value 
of 20.4 K for a charge of 0.5 electron per unit cell, which 
compare well with the experimental value. 



sponse theory, for the electron-phonon coupling constant 
A, the weighted average of phonon frequencies oJiog^ and 
the density of states at the Fermi level Np. We have used 
the McMillan formula to estimate Tj^ 



T n = 



E20 



exp 



-1.04(1 + A) 
A-^*(1 + 0.62A) 



(1) 



with n* = 0.1 to obtain an agreement between our cal- 
culations and the reported T c of 11.2 K for NbC Pit is 
important to note that this equation is valid for A < 1. 
In this context, the reported T c for NbN 5 is inaccurate 
as A is too large for the McMillan approximation. Note 
that the McMillan formula can be used with different 
definition of the prefactor (w) . We used the definition of 
w iog = exp (In w) because it gives a better agreement be- 



IV. DISCUSSION 



The coupling constant A can be separated into two 
contributions, one from the optical modes and the other 
from the acoustic modes denoted A ac . The values of A ac 
are reported in Table [T] together with A. We can see that 
the variation of A as the density of states increases comes 
mainly from the change of A ac while the optical part con- 
tributes only a constant independent of Np . This is con- 
sistent with the hypothesis that the enhancement of the 
Kohn anomaly dominates the evolution in the electron- 
phonon coupling. Following this idea, we propose a sim- 
ple model to express A and T c as functions of Np. 

The values of A and wi og can be extracted from the 
Eliashberg spectral function a 2 F(uj)^ This function is 
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defined 

a 2 F(u) = ^ ]T a 2 F ( k > k ', w ) S " zf) S (e' fc - e F ) 

F kk' 

(2) 

a 2 F(k, k', u) = J2 N FW k ,k> | 2 <5 (w - Wk-k'j) 

<4, = (h/2M^u^) 1/2 (fc|e k _ k ,j • W|fc') 

where e qj - is the normalized atomic displacement for a 
phonon. The sum over j is over the different phonon 
branches. Since g is inversely proportional to y^cj, the 
amplitude of the movement of the atoms is enhanced near 
the structural phase transition where u becomes small for 
a soft mode, resulting in a stronger coupling. With all 
this in mind, we can represent the a 2 F(u>) function by 
setting apart the contribution of the softened mode at 
frequency u': 

a 2 F(u) = N F f (cj) + m2 ^ F S(uj - J) , (3) 

U)' 

where \M\ 2 represents an effective coupling matrix ele- 
ment of the relevant soft modes. From Green's function 
theory of electron-phonon coupling,^ the renormalized 
phonon frequencies can be expressed in terms of the elec- 
tron polarization, n(w, q) as 



{l + Tl(LO,q)g 2 



(4) 



where g is the coupling and loq is the bare phonon fre- 
quency. Since the electron polarization is proportional 
to —Np, the renormalized phonon frequency of the form 
uj' = loq yjl — Nf/N c is used to describe the evolution of 
the softened phonon, where N c is the critical density of 
states at which the structure becomes unstable. 

The function fo(to) is independent of Np and repre- 
sents the contribution from all the phonons that are not 
affected by the Kohn anomaly. Using ([3]), we obtain an 
expression for A: 



A = 2 / duj 



a 2 F{uj) 



= V N F 



2\M\ 2 N F 



(5) 



where Vq = 2 J du>fo(uj)ui 1 is a material dependent pa- 
rameter. This also gives a form for wi og 



exp 



CN F 



\M\ 2 N F 



In a/ 



(6) 



where C — J dicfoiuS) ln(aj)a; _1 is another parameter. 
Although this is a simple model for the Eliashberg spec- 
tral function, it allows a qualitative description of the 
effect of the soft phonon on the T c . A possible improve- 
ment of the model would be to use a localized function in 
the second term of Q to represent the phonon softening. 
For example, we could use a gaussian function to repre- 
sent the softening which will allow to better describe the 



1 1 1 1 1 1 1 1 1 1 1 1 

□ Computed values (NbC [ 


ii |. i 


A Computed values (NbN + charge) 




- — Model (NbC ] x N x ) 




■ ■ ■ Model (NbN + charge) 


A 


□ 




a ' 


A 






1 ... 1 ... 1 


1 



0.38 0.4 
N p (states eV 1 cell ') 



FIG. 4. (color online). Fit solutions for X(Nf) based on 
|5| for mixed pseudopotential method (black curve) and for a 
uniform background charge (red curve). Parameters are fitted 
to the computed values and are reported in Table |H| 



phonons affected. In order to simplify the model, a Delta 
function is used. 

We can express X(Np) by using only three parameters, 
\M\ 2 ojq 2 , N c and Vq. Our calculations of the phonon 
spectra yield an approximate value of N c . We used a fit 
for the other parameters which are reported in Table |H| 
Results for these parameters for the case of NbCi-^N^ 
and the nitrogen deficient NbN are similar. 



TABLE II. Fit solution for the parameters in |5| and j6j. 
Contribution of the normal phonons to the electron-phonon 
coupling are represented by Vq. The softened frequency is 
taken into account by \M\ 2 which is the coupling matrix ele- 
ment and by wo, the bare frequency for these phonons. The 
critical density of states, N c , indicates the occurrence of the 
structural phase transition. The parameter C is obtained by 
inverting |6| and by using the values reported in Table [T] 





NbCi-sN* 


NbN 1 


2\M\ 2 lJq' 2 (eV cell) 


0.023 


0.025 


iV c (states eV" 1 celP 1 ) 


0.414 


0.424 


V (eV cell) 


1.87 


1.41 


C (eV cell log(eV)) 


-3.38 


-2.43 



The model for \{Np) and the computed values are 
shown in Fig. |4j To compute T c using the McMillan 
formula (JTJ , the bare phonon frequency is needed in order 
to determine LU\ og . Based on Fig. [3j the renormalized 
frequency for the soft phonon is estimated at 20 meV for 
NbC. Hence, the bare frequency is estimated by luq — 
w (1 — Np/N c ) 1 m 60 meV. This gives an expression 
for T c (Np) which is plotted along with the calculated 
values in Fig. [3J Our model predicts an increase in T c as 
Np increases and a sharp fall when Np rj N c . However, 
this feature appears in the region where A is very large. 
It is known that the McMillan equation is not adequate 
if A > 2. In a recent work lower and upper bounds for 
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tween the weak-coupling limit for T c cx exp(— 1/A) and 
the strong-coupling one T c cx 

When used in our model, these bounds do not fall to 
zero at Np ~ N c where A —> oo. Using the correct limit 
of T c for large M. 



T c 



'Xuj. 



ph 



(9) 



where ui p h is the average of phonon frequencies which can 
be taken to be uj\ g- Near the phase transition, only the 
contribution of the soft frequency is important. Hence, if 
Np ~ N c , according to ([5]) and ([6]), we find that A ~ uj'~ 2 
and cjjog ~ ui' so that {f9j) gives T c ~ y/2\M\ 2 N c which is 
indeed not zero. 



FIG. 5. (color online). Best fit solutions for T c (Nf) based on 
the McMillan equation Q. The cji og term was computed us- 
ing pL fi* was fixed at 0.1 in order to get agreement between 
our results and the experimental T c for NbC™T c falls rapidly 
to zero near the critical DOS which is the onset of the struc- 
tural phase transition. This is due to the inadequacy of the 
McMillan formula when A > 2. The experimental values of 
T c reported on this figure are taken frorrP assuming that the 
relation between nitrogen concentration and the calculated 
DOS with alchemical pseudopotential is adequate. 



T c have been derived: 

/(A(c 2 )/ H 2 nax ,M*) ^ max < T c < f (A, fj,*) (7) 
where 

is an empirically smooth function that interpolates be- 



V. CONCLUSION 

In conclusion, we have presented ab initio pseudopo- 
tential DFT results on NbCi-ajN^ based on the vir- 
tual crystal approximation. We showed that the Kohn 
anomaly leads to an enhanced electron-phonon coupling. 
These results can be modeled by singling out the soften- 
ing contribution to the Eliashberg spectral function. 
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